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Abstract 


A systematic and comprehensive investigation was performed to provide detailed data on the 
three dimensional viscous flow phenomena downstream of a modem turbine rotor and to understand 
the flow physics such as origin, nature, development of wakes, secondary flow, and leakage flow. The 
experiment was carried out in the Axial Flow Turbine Research Facility (AFTRF) at Penn State, with 
velocity measurements taken with a 3-D LDV System. Two radial traverses at 1% and 10% of chord 
downstream of the rotor have been performed to identify the three-dimensional flow features at the exit 
of the rotor blade row. Sufficient spatial resolution was maintained to resolve blade wake, secondary 
flow, and tip leakage flow. The wake deficit is found to be substantial, especially at 1% of chord 
downstream of the rotor. At this location, negative axial velocity occurs near the tip, suggesting flow 
separation in the tip clearance region. Turbulence intensities peak in the wake region, and cross- 
correlations are mainly associated with the velocity gradient of the wake deficit. The radial velocities, 
both in the wake and in the endwall region, are found to be substantial Two counter-rotating 
secondary flows are identified in the blade passage, with one occupying the half span close to the 
casing and the other occupying the half span close to the hub. The tip leakage flow is well restricted to 
10% immersion from the blade tip. There are strong vorticity distributions associated with these 
secondary flows and tip leakage flow. The passage averaged data are in good agreement with design 
values. 
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Nomenclature 


Cr rotor blade axial chord at mid-span 

H normalized distance between the hub and the casing. (r-r h )/(r t -rh) 
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mass flow rate 

number of measurements per bin 
pressure surface 
radius of the hub 
radius of the blade tip 

streamwise, normal, and tangential directions 
tangential distance normalized by blade spacing 
suction surface 

the blade mid-span rotational velocity 

fluctuating velocity components (normalized by local velocity) 

absolute velocity 

velocity measurements in the LDV system coordinate directions (blue, green, and 

violet channels). 

relative velocity 

fluctuating relative velocity 

shear stress 

wake deficit velocity (w^-w,™) 

axial, radial, and tangential direction (x/Cr = 0 at LE; = 1 at TE) 

normalized axial distance (x/Cr) 

flow angle measured from the axial direction 

radial flow angle measured from axial direction 

LDV channel angles 

axial vorticity (Eq. 3) 

the rotor angular velocity 


Superscripts and Subscripts 

e ensemble average 

BA blade aperiodic component 

BP blade periodic component 

max maximum 

min minimum 

RP revolution periodic 

RA revolution aperiodic 

s,n,r streamwise, normal, and radial directions 

x,r,0 axial, radial, and radial direction 

1 ,2,3 upstream of nozzle, upstream of rotor, downstream of rotor positions 
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1. INTRODUCTION 


The flow field in an axial-flow turbine rotor is highly three dimensional. The distinctive three 
dimensional viscous effects include blade and endwall surface ( laminar , transitional and turbulent) 
boundary layers, shock-boundary layer interaction, blade wake, horseshoe vortex near the blade leading 
edge, tip leakage flow, and the secondary flow. 

Much of the research work on turbine rotor flows has been confined to measurements in 
cascades, and substantial understanding of the three-dimensional flow field in the blade passage has 
been gained. The pioneering work by Langston et aL (1977) of secondary flow patterns in a turbine 
cascade gave a detailed description of the passage vortex and the leading edge horseshoe vortex. 
Sieverding (1985) summarized extensive results of experimental investigations on the synchronous 
evolution of horseshoe and passage vortices and on the effect of vortex structures on end-wall 
boundary layers. The consensus was that the pressure side branch of the horseshoe vortex followed a 
smooth curve through the blade passage as part of the passage vortex, while the suction side branch 
wrapped itself around the passage vortex. Gregory-Smith (1988) and his group have acquired valuable 
data in the secondary flow region in a cascade, including turbulent stresses, intensities, and losses. 
Concerned with the mixing loss downstream of a turbine cascade, Moore et aL (1985) found that the 
rise in total pressure loss was almost entirety balanced by a decrease in the secondary flow kinetic 
energy. The tip leakage flow was, on its own complexity, the subject of many turbine cascade 
researches, and a recent work by Yamamoto (1989) in a cascade demonstrated the strong effect of gap 
size and inlet incidence on the leakage flow and the passage vortex. Focused on the convective 
transport from a simulated turbine blade in a linear cascade, Goldstein et aL (1995) found that the 
endwall flows enhanced the naphthalene sublimation substantially and that the disturbance near the 
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leading edge would trigger an early laminar-turbulent transition. 

Due to the general industry trend towards highly three-dimensional blading geometry to rreet 
the requirement of increased aerodynamic efficiency with fewer tuibomachine components, 
investigations related to flow conditions in real turbine stages have been intensified. Joslyn et aL (1982) 
examined the unsteady three-dimensional flow near the mid-span of the blade in a one-and-half stag e 
axial turbine and found that, downstream of the rotor, the pitch angle of the flow even exceeded the 
acceptance angle of their three-element hot-wire probe during part of the blade period. Later, Joslyn et 
aL (1992) measured the mean properties of the three-dimensional flow field in the saire turbine 
configuration using the five-hole probe area traversing in the rotating flame. Their measurements were 
further complemented by the surface flow visualization. They found that the hub and tip secondary 
flows, as well as the relative eddy and the tip leakage flow, were the main causes for the three- 
dimensionality of the turbine rotor flow field and that the pressure distribution on the rotor blade was 
strongly influenced by this three-dimensionality. Hodson and Addison (1989) and Hodson (1993) 
provided comprehensive steady and unsteady flow data at the micspan of a low-speed rotor, including 
boundary layer and loss characteristics and transition. Gallus et aL (1995) pointed out that the 
convection of stator secondary vortices into the rotor flow field had received much less attention than 
the other two forms of rotor-stator interaction (the potential flow field and the blade wake). Their 
experimental and numerical results demonstrated a much improved prediction of the rotor passage 
vortex and the flow exit angle near the hub by considering the unsteady interaction between the stator 
secondary flow and the rotor flow field. They also found that the tip leakage vortex is confined to a 
region very close to the shroud due to a very tight tip clearance. 

The objective of this study is to provide comprehensive data about the three dirransinnai 
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viscous flow phenomena downstream of a modem turbine rotor blade passage and to understand their 
physics such as origin, nature, interaction, and effect on the overall performance. The velocity 
measurements have been carried out in the turbine facility with a 3-D Laser Doppler Velocimeter. 
Two radial traverses downstream of the rotor have been recently accomplished to identify the three 
dimensional flow phenomena in the wake of the rotor blade. 

2. FACILITY AND INSTRUMENTATION 

The Axial Flow Turbine Research Facility (AFTRF) at Penn State is an open circuit, single 
stage with a diameter of 91.4 cm and a hub-to-tip ratio of 0.73. It consists of a bellmouth inlet, a 
turbulence generating grid section, and a test section consisting of a row of 23 nozzle guide vanes and 
a row of 29 rotor blades. The configurations of the vane and blade are designed by the GE Aircraft 
designers and represent the features of a current or near-future HP turbine stage (Halliwell, 1991). 
Following the test section are two 100 hp, adjustable pitch, axial flow fens in series to provide a 
through flow up to 22,000 cfin. The power generated by the turbine rotor is absorbed by an eddy- 
current brake, and the rotational speed of the rotor can be varied between 175 and 1695 rpm. The 
peak total-to-total isentropic efficiency is 89.3% at a rotational speed of 1300 rpm. At this peak 
efficiency point, the total pressure and temperature ratios are P03 /P01 = 0.928 and T 03 /T 0 i = 0.981. 
(For more details, see Lakshminarayana et aL, 1992.). The experiment was conducted at a corrected 
rotational speed of 1300 rpm and a mass flow rate (m ) of 10.53 kg/sec. Design features of the 
AFTRF are shown in Table 1. 

The T as&r Doppler Velocimeter (LDV) system is a three-color cross-beam system (TSI Model 9100- 
12). The li ght source is a 7 W argon-ion Laser (Coherent Model Innova.300), operating in the multi 
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line mode. As shown in figure 1, the blue and green light beams are introduced into the flow f Hd. in 
one optical rail, perpendicular to the glass window on the turbine casing to sample axial and tangential 
components of the velocity. The violet beam is arranged in an another optical rail with an ang le of 
27.8° from the perpendicular in order to measure the radial component. In addition, 32.5° and 45° 
angle arrangements are also tried in the measurements, and it is found that the configuration with 
around 30° angle gives a sufficient resolution of the r adial velocity while still keeping the system 
manageable. The displacement caused by passing a light beam at an incidence angle through a glass 
window is accounted for by placing a similar glass in the optical path when the beam crossing 
alignment is being made. Forty MHz frequency shifting are used to reduce the fringe bias and to solve 
the direction ambiguity by passing one of the two crossing beams through a Bragg celL 
Thus, the three dimensional velocity can be calculated from: 

V, = Vcsiny, -i- VySinYi 
an(y, + y 2 ) 

„ _ Vo cosy, - VveosY, 

V r — ; - " vlj 

sin(Y,+Y 2 ) 

V e = V B 

where V B , Vo, and Vv are the velocities sampled by the blue, green and violet light beams, respectively 
and Yi and Y 2 are the angles of the green and violet optical paths. In the present investigation Yi = 27.8° 
and Yi = 0°. 

The receiving optics are arranged off-axis, i.e. the photo multiplier for the violet beam receives 
scattered light from the seeding particles in the blue and green trans mitting optical r ail, and vice versa. 
This arrangement substantially reduces noise from internal refection and li ght reflected and scattered 
from the glass window. Liquid particles with a nominal mean particle si7« of about 0.6 fim, generated 
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Table 1 The Design Features of AFTRF 


Hub Tip Ratio 

0.7269 

Tip Radius 

0.4582 m 

Blade Height 

0.1229 m 

Tip Rel. Mach Number 

0.240 (max) 

Rotor Blade (tip) 


number 

29 

chord 

0.1287 m 

spacing 

0.1028 m 

turning angle 

95.42° (tip), 125.69° (hub) 

maximum thickness 

22 mm 

tip clearance 

1.27 mm (actual 0.97 mm av, 1.04 mm max, 0.77 mm min) 

Nozzle Guide Vane (tip) 


number 

23 

chord 

0.1768 m 

spacing 

0.1308 m 

turning angle 

70 deg. 

maximum thickness 

38.81 mm 

Auxiliary Fan 


pressure rise 

74.72 mm Hg 

mass flow 

10.39 m 3 per sec. (22,000 cfin) 

power 

149.1 kw 

Vane Reynolds Number 


based on inlet vel 

(3 -4) x 10 s 

based on exit veL 

(9 - 10) xlO 5 

Blade Reynolds Number 


based on inlet reL veL 

(2.5 - 4.5) x 10 s 

based on exit reL veL 

(5 - 7) xlO 5 
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by a 6-jet atomizer (TSI Model 9306), are introduced into the flow field at the bellmouth inle t of the 
facility as the seeding particles. In selecting the working liquid for the atomizer, propylene glycol is 
decided on largely due to its low toxicity and water solubility. The latter is especially helpful for the 
investigations of the tip leakage flow, since the seeding particles directly contaminate the glass w indo w 
and the water solubility allows for a quick clean-up during the experiments. The whole LDV system is 
mounted on an optical table, and a hydraulic system is provided for the table movement in x, y, and z 
directions, plus a tilting. Thus, traversing of the LDV probe volume in radial, tangential and ayial 
directions can be readily achieved. 


3. DATA ACQUISITION AND DATA REDUCTION PROCEDURE 

The current LDV set-up, with the violet beam inclined at 27.8° from the perpendicular on the 
downstream side of the blue and green beams (see Fig. 1), is arranged ideally to acquire data 
downstream of the rotor for the least shadow effect from the blade. A PC-based interface (Dostek 
Model 1400 A) records the three velocity component data along with the rotor orientation from a shaft 
encoder (6000 pulses per revolution clocking and 1 pulse per revolution tri ggeri ng) All m easurements 
were acquired with the facility operating at peak efficiency. 

The data reduction is accomplished by ensemble averaging and then by passage-to-passage 
averaging according to the recorded rotor orientation. The entire rotor is discretized into 1450 bins 
(50 per blade passage as can be seen in Figure 2.) and the corresponding measurements are averaged 
over the individual bins to yield the ensemble average velocity. One of the main limitations of the LDV 
technique is apparent in the non uniform distribution of measurements across the individual blade 


6 



passages shown in Figure 3 and 6. A fairly low seed density in the secondary flow region required an 
enormous amount of measurements to resolve the blade-to-blade variations accurately. To determine 
the minimum data number of data points per window required for a sufficient passage averaging, we 
examined the mean velocity and the turbulence intensity (the rms value of the instantaneous velocity 
minus the passage-averaged mean) at various stations with sample sizes from 50,000 to 200,000 as 
shown in figure 4. The passage-averaged properties were computed and compared for a variety of 
individual bin measurements and found to be essentially constant for samples larger than 100 
measurements per bin. 

In addition, the actual number of data windows per blade passage, each representing a finite 
time interval, was varied to determine the optimal size. If within this time interval, there exists 
appreciable variation in velocity, especially in the wake region, artificial fluctuations may influence the 
averaging procedure and result. The number of bins per blade passage was varied between 50 and 100 
(as shown in fig. 5) at a few selected stations with large velocity gradients and negligible variations 
were seen between the two temporal discretizations. A sample size of 50 was chosen to keep the 
overall data sample size manageable, while at the same time ensuring measurement accuracy. 

The following data are derived, presented, and interpreted in this paper. 

1 . Three components of mean velocity is cylindrical co-ordinate system, (r, 0, z) 

2. Secondary velocity vectors and axial vorticity 

3. Total turbulence intensity and stresses in the intrinsic co-ordinate system (s, n, r) 

The intrinsic co-ordinate system (s, n, r) is also used to derive the mean velocities, turbulence 
intensities, and stresses. This provides valuable information on the stresses in the streamwise and radial 
directions. In this co-ordinate system, s is the streamwise direction, n is the principal normal (blade-to- 
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blade) direction and r is the radial direction. The s co-ordinate is aligned with the de sign flow angle 
tabulated in Table 2. Any deviation from this angle is the secondary flow represented by 

Vsec = V design “ V measircd (2) 


Table 2 The design value of the relative flow exit angle at the blade trailing edge. Here, r t is the 
radius of the blade tip, and n, the radius of the hub. 


H 

1.0 

0.9 

0.8 

0.7 

0.6 

0.5 

0.4 

0.3 

0.2 

0.1 

0.0 

-P(°) 

68.15 

67.87 

67.55 

67.26 

66.94 

66.64 

66.32 

65.97 

65.64 

65.32 

65.00 


The mean velocities presented in this paper are based on the “averaged passage” derived from 
averaging the ensemble averaged velocity across all the blade passages. This procedure is considered 
valid for the investigation of single-stage flow field of the AFTRF. To verify that the aver aging 
procedure is valid, the procedure developed by Suder et aL (1987) and later modified by Suryavamshi 
et aL (1994) to decompose the instantaneous velocity measurement in a multi-stage environment was 
used at a few selected measurement stations. The decomposition is based on the following: 

V = V + V BP + Vba + Vra + v' (3) 

where: V is the time mean, Vbp is the blade periodic, V B a is blade aperiodic, Vra is the revolution 
aperiodic, and v' is the unresolved unsteadiness due to turbulence and revolution-to-revolution 
variations. The summation of V BP and V B a represents the revolution periodic (Vrp) component. The 
general procedure in the computation of the individual components involves computation of the time 
mean ( V ) and the ensemble average (V e ) over the entire turbine rotor. The revolution aperiodic (Vra) 
component is computed as the average variation of the difference between the time mean and ensemble 
averaged velocity over the individual blade passage. For a single-stage facility the Vra component is 
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essentially zero. The revolution periodic (Vrp) component is computed by subtracting the Vra 
component from V e - V at each bin location. The Vrp component can now be decomposed into the 
blade periodic (V B p) and blade aperiodic (V B a) components. The blade periodic is representative of the 
flow as seen by the “average” blade passage and the blade aperiodic component is a measure of the 
blade-to-blade passage variation. Figure 3 shows the relative magnitudes for the axial velocity 
components in the first six blade passages at x/Cr = 1 .01 and H = 0.93 locations. As expected, the Vra 
velocity component is negligible (maximum of 0.02U,,,). The blade aperiodic component is relatively 
small over most of the blade passage except in the wake region. 

To assess the fluctuation in velocity caused by the variation between individual blades involved 
in the passage averaging process used in this study, a large amount of data is sampled at one location 
(x/Cr = 1.10 and 95% immersion) with at least 1500 data points per window for the average passage in 
the wake is utilized. Thus, the minimum data number per window for each individual blade passage is 
around 1500/29 = 51. As a demonstration of this argument, the axial turbulence intensity for five 
individual blade passages (based on the ensemble-averaged mean) and for the average passage are 
plotted in Fig. 7. The clustering of individual blade passages data around the average passage suggest 
that the turbulence intensity distribution is nearly identical in all passages. 

This clearly demonstrates that the passage averaging procedure can be used to derive the 
blade-to-blade variation of the three components of velocity. The passage-average properties are then 
used to derive the unresolved unsteadiness in the rotor passage using Eq. 3. The unresolved 
unsteadiness at x/Cr =1.10 and H = 0.6 is shown superposed on the blade periodic quantities in Fig. 8. 

Two radial traverses of LDV measurements (a total of 18 radial locations indicated in Table 3, 
with higher resolution near the blade tip) were completed at 1% and 10% chord length downstream of 
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the rotor trailing edge (denoted as x/Cr= 1.01 andx/c T = 1.10) to identify the various flow features. 

Table 3. Measurement locations in the radial direction at x/Cr=1.01 and 1.10. 


H = 

0.02 

0.05 

0.10 

0.20 

0.30 

0.40 

0.50 

0.60 

0.70 

0.75 

0.80 

0.85 

0.90 

0.93 

0.95 

0.97 

0.98 

0.99 


4. RESULTS AND DISCUSSION 
4.1 Blade-to-Blade Profiles 

The mean axial, radial, relative tangential, total velocity, and relative flow angle are shown in 
Figures 9 through 13 for all the radial measurement locations at the x/c, = 1.01 station. S imilar hlade to 
blade profiles for the x/c,=1.10 measurement location are shown in Figures 14 through 18. Turbulence 
intensities in the s,tur coordinate system and the respective correlations are shown in Figures 19 
through 30 for the two downstream locations. 

The mean velocity profiles at x/c, = 1.01, shown in Figs. 9 through 13, clearly indicate the 
presence of thick boundary and wake on the suction side and thin boundary layer and wake on the 
pressure side. The fact that the suction surface boundary layer wake is substantially thicker than those 
on the pressure side indicate that the probable cause of this is the transport of nozzle wake towards the 
suction side (unlike compressor). This results in higher velocity defect and higher turbulence 
intensities. Near the casing (H= 0.93 and 0.98) the flow is nearly hilly developed, with substantial r adial 
inward velocity near the suction surface caused by the secondary flow. Very strong secondary flows 
and leakage flows interact in this region. The flow in the secondary flow regions is overturned and the 
turbulence intensities are high. The near wake features, where the turbulence intensities are lowest in 
the wake center, can also be clearly seen. Furthermore, the radial and normal components (Figure 19- 
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21) of turbulence intensities are higher than the streamwise component in the wake and secondary flow 
regions, thus confirming the effect of rotation in modifying the turbulence structure. The Coriolis 
forces tend to increase the turbulence intensities in the radial direction inside the rotor boundary layer, 
and this is consistent with the turbulence measurement and analysis inside an inducer (Anand & 
Lakshminarayana, 1978) and at the exit of a compressor rotor (Lakshminarayana and Reynolds, 1979). 
The shear stresses in the streamwise direction (w s - w n 0 shown in Fig. 22 are higher than those in the 
radial direction. 

The data away from the casing (H = 0.5 in Fig. 9 through 13) shows that the suction surface 
boundary/wake thickness decreases towards midspan, being fairly thin with low radial velocities at the 
mid-radius. One of the interesting features observed at H = 0.8 is the presence of thick wakes and high 
turbulence intensities. The probable cause of this is the presence of nozzle secondary flow, which was 
observed at this location (Zaccaria and Lakshminarayana, 1995). The data in the hub region (at H = 0.1 
in Fig. 10) indicates the presence of very large radial outward velocity, with a thick wake on the 
suction side comer. The wake description is incomplete due to a lack of accessibility in this region by 
the Laser Doppler Velocimeter. 

The freestream turbulence intensities, shown in Figures 19-21 and 25-27, vary from 4 to 6% 
across the span, with the exception of locations very close to the hubwall region. The turbulence 
intensity at inlet to nozzle vanes varies from 1 to 1.5%. Hence the free stream turbulence has 
increased substantially due to the interaction of the nozzle wake with the rotor flow field. 

The region from 20 to 60% span (Figs. 9-18) is found to be dominated by inviscid core (where 
the flow is uniform) and the blade wake. The thinnest wake occurs at H = 0.5 and the wake becomes 
increasingly wider toward the endwall regions. No core or inviscid flow region is found near the 
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casing or hubwall regions. 


4.2 Secondary and Leakage Flow Field 

The secondary velocity vectors, defined by Eq. 2, are shown in Fig. 31. Five distinct regions of 
the flow field can be identified; passage core flow, hub passage vortex, casing passage vortex, rotor 
wake flow, and the tip vortex flow. The blank region near the hub blade region at x/c, = 1.10 is due to 
inaccessibility of the flow field by the laser. It is clear that there are two secondary flow regions in the 
blade passage. Toward the hub, a strong, well-organized and smooth secondary flow pattern 
occupying the entire half-span is observed, due to large blade turning angle. Its center is at 
approximately 20% pitch from the blade suction side and between H = 0.30 and H = 0.20 immersion. 
However, the secondary flow and leakage flow in the casing half of the passage is much more 
complicated and much stronger. The radially outward flow on the blade pressure side possesses 
significant radial outward velocity only when it is very close to the casing (within 13% of the blade 
passage). The radially inward flow on the blade suction side, nevertheless, covers the whole half-span. 
The core of the secondary vortex is located near H = 0.85 and S = 0.4. The secondary flow tends to 
hug the tip leakage vortex and contains it in the suction surface - casing comer. The leakage flow 
phenomena observed here is substantially different from those observed in cascades (e.g. Yamamoto, 
1989). Unlike cascades, the tip clearance vortex is confined to a small region near the suction surface - 
casing comer. In cascades, the leakage flow has a tendency to move much further away from the 
casing and away from the suction surface. The secondary flow and vortex has a much stronger 
influence than the leakage vortex for this particular rotor. 

Because of the leakage flow, the secondary flow pattern in the casing-half of the passage 
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becomes much distorted. The leakage and secondary vortices entrain fluid from surroundings as 
confirmed by outward velocity near the pressure side and inward velocity on the suction side. The 
secondary flow is confined to the region above H = 0.90. Below H= 0.90, high radial velocity (inward) 
occurs in the wake - secondary flow interaction region. This region has radial inward velocity and 
involves interaction and mixing of rotor wake and secondary flow. It is characterized by the radial 
flows, low total velocity and high vorticity as explained later. Thus, it is clear that rotor wakes are 
substantially influenced by secondary flow from mid-radius to tip, with secondary flow effects 
dominating (over rotation effects) in inducing inward flows. The wake region from hub to mid-span 
region has radial outward velocity, this is consistent with the effects of rotation on blade boundary layer 
and wake. The rotation and secondary flow effects have a cumulative effect in inducing large radial 
outward flow in the wake regions from the hub to mid-radius. 

Further downstream at x/Cr = 1.10 (see Fig. 31b), the leakage flow weakens and is nearly 
circular. The secondary flow becomes dominant near the blade suction side-casing comer. The zero 
r adial flow velocity region on the blade suction side moves inward up to H = 0.85. The rest of the flow 
features remains the same but at a slightly weaker level. One of the interesting features of this flow 
field is that the growth and diffusion of leakage vortex are substantially less than those observed in 
cascades. The vortex tends to stay near the suction surface comer closer to the casing and its extent is 
reduced. This is brought about by strong secondary flows, which counteract and prevent growth and 
spreading of the leakage vortex. 

43 Contours of Mean Velocity Components 

The contour plots of axial, relative tangential velocity, relative flow angles and radial angles at 
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x/Cr = 1.01 and 1.10 are plotted in Figs. 32, 33, 34, and 35 respectively. The blade-to-blade profiles are 
shown in Figs. 9-18. In the midspan region, the wakes are thin and axial velocities are high. The 
secondary and leakage flow regions have the complex flow characteristics. The wake width increases 
continuously from mid-span to hub and mid-span to casing; the endwall region (H < 0.1, H > 0.9) 
shows nearly fully developed flow. The tip clearance region has separated flow, with substantial 
negative axial velocities. The secondary flow region (around H = 0.85) shows reduced axial velocities. 

The blockage effects due to secondary and leakage vortices/flows are also evident from this plot. 
The region away from these flows shows very high axial velocities extending all the way from H = 0.6 
to H = 0.9. As the flow progresses downstream from x/Cr = 1.01 to 1.1, the reversed flow disappears, 
the wake and secondary flow region widens, and the highest axial velocity occurs near the midspan. 
But for these changes, the distributions is s imilar to that observed at x/Cr = 1.01. 

The relative tangential velocity (WJ is shown in Fig. 33. The tip vortex has low axial and 
relative tangential velocities. In the adjoining secondary flow region, the relative t angential velocity 
increases rapidly away from the blade, indicating the usual overturning caused by secondary flow and 
vortex. The core of this vortex is located near H = 0.85 closer to the suction side. Below this core the 
relative tangential velocity decreases, indicating flow undertuming in this region. This confirms the 
earlier conclusion that the secondary flow dominates the outer span of the rotor. 

The design value of the blade exit angles are listed in T able 2 and the measured values are 
shown in Fig. 34. As demonstrated in the contour map, in the core inviscid flow region, the measured 
flow angle is very close to the design value. In the wake region, a distinctive pattern exists. On the 
blade suction side, the flow is overturned, and on the pressure side, the flow is undertumed. Near the 
mid-radius the relative flow angle changes from -79° on the suction side to -56° on the pressure side. 
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Near the hub, especially at the blade suction side-hub comer, the flow is overturned as a result of both 
the strong secondary flow toward the blade suction side (see Figs. 31 and 34). The flow overturning 
caused by the casing secondary flow at the blade suction side-casing comer is also visible, but 
overpowered by the tip leakage flow. The strong leakage vortex and associated very low axial velocity 
results in the flow angle approaching -90°. In feet, at x/c, = 1.01, negative axial velocity exists in the tip 
leakage vortex and the flow angle is beyond -90°. The influence of tip leakage and secondary flows on 
the flow angle is so strong and so extensive that the previously described wake pattern of the flow 
angle distribution (overturning on the blade suction side, undertuming on the pressure side) is 
completely distorted for H > 0.7 (see Fig. 34). On the blade pressure side, the tip leakage flow causes 
flow undertuming in a very confined region near the tip gap. 

Near the midspan region (H = 0.4 through 0.85) the axial and radial velocity variation is very 
small outside the wake region. The W e variation across the blade passage is nearly linear, influenced by 
the blade-to-blade pressure gradient. In the wake, the velocity defect is quite large accounting for a 
drop in the axial velocity from a free stream value of 0.7 U m to a minimum of only 0.18 U m in the 
wake. Toward the tip region, the wake width increases due to the casing suction passage vortex 
interacting with the low momentum fluid in the wake region. At the same time, the radial velocity 
variation grows dramatically from almost negligible levels in the midspan locations to large radially 
inward levels of -0.42 U m at H = 0.7 (Fig. 31). Toward the tip, the radial velocity magnitude starts to 
drop, as the secondary flow influence decays, till about the H = 0.93 location. 

In the tip clearance region extending from approximately H = 0.95 to the casing, the radial and 
axial velocity profiles start to exhibit very large gradients due to the complex interaction between the 
casing endwall boundary layer, tip clearance vortex and the blade wake. The tip clearance vortex is 
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clearly visible in the radial velocity profiles at H = 0.95 and H = 0.97 with very high velocity gradients 
on either side of the vortex. The tangential velocity profiles in this region still exhibit a large li near 
gradient region in the free stream and a relatively small wake region when compared to the axial and 
tangential velocity profiles. 

The radial flow angle, shown in Fig. 35, remains close to zero in the mid-portion of the 
blade (H=0.60 and H=0.50). Toward the hub, encompassing the hub passage vortex region a clear 
distinction can be made between the radially outward flow on the blade suction side and inflow on 
the pressure side. In the casing vortex region, the trend is opposite, with an inflow on the blade 
suction side and an outflow on the pressure side. Very close to the casing (H=0.99), however, the 
radial velocity drops close to zero again. 

The measured velocity and angle distributions at 10% chord downstream exhibited similar 
behavior to that at the 1 % location (Figs 3 1 - 35). The major difference is in the suction-casing corner, 
where the overall level of axial velocity decreased toward zero compared to the negative velocity 
region at x/Cr=1.01. The decay is due to mixing and diffusion caused by the radial outward motion of 
the tip clearance vortex and an increase in the wake width in the downstream direction. The overall 
flow features remain very similar, indicating a strong and well established secondary flow pattern 
propagating downstream In general, the following conclusions can be drawn with regard to the flow 
variation in the downstream direction x/Cr= 1.01 to 1.10: 

1. The radial velocities are slightly lower at the downstream location due to miring and 
decay. 

2. The hub passage vortex remains well organized and at the same position relative to the 
blade passage. 
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3. 


The casing passage vortex moves toward the suction side of the wake due to the low 
pressure region created by the tip clearance vortex and the overall blade-to-blade 
pressure gradient. There is very little movement of the passage vortex in the r adial 
direction as the flow progresses downstream. 

4. The tip-clearance vortex remains well organized, decaying in strength and moving 
slightly inward radially due to r adial pressure gradient and the large scale casing 
passage vortex motion. 

The radial distribution of the passage-averaged velocity components at x/Cr=1.10 is compared 
with the design values in Fig. 36. The averaged radial velocity remains very close to zero, even though 
there are strong secondary and tip leakage flows existing in the blade passage. The measured axial and 
total velocities are lower than the design value in the hub-half of the span due to stronger secondary 

flow there. In the outer span region, the total velocities are slightly higher. The casing region (90 - 

100%) has lower absolute tangential velocity (over turning). 

4.4 Axial Vorticity 

The axial component of vorticity can be derived exactly since all the velocity gradients can be 
calculated in (r, 6) directions from the measured data. The gradients in the axial direction are not 
accurate due to the coarse measurement grid. The axial vorticity is calculated from the equation 

©* = 1 9(rWe) _ iWr (3) 

£2 r dr rd0 

The contours of axial vorticity are shown in Fig. 37. The negative values indicate anti clockwise 
motion. The leakage vorticity is confined to the casing within 5% of the span. The wake vorticity is 
positive (due to inward radial velocity) in the outer span. The secondary vorticity (negative) near the 
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casing region is a dominant feature of the flow field as indicated earlier, and it covers substantial 
portions of the passage. Likewise the secondary vorticity is very strong near the hub region and the 
m axi m u m strength is between the hub and midspan region, closer to the suction surface. The strength 
of the maximum leakage vorticity, secondary vorticity in the hub and endwall regions are about the 
same orders of magnitude. Furthermore, the wake-secondary flow interaction region near 60 -70% 
span, close to suction surface, indicates the presence of substantial vorticity. 

The highest level of this vorticity distribution occurs at H = 0.35 immersion and 15% pitch 
from the blade suction side, which is also the region of high turbulence intensity (as explained later). 
For the casing secondary flow, however, a negative vorticity distribution exists in a more fragmented 
form due to the interaction between the secondary flow and the leakage flow. On the suction side and 
within 15% radial immersion, the vorticity is negative. The positive vorticity associated with the 
leakage vortex is clearly seen at the blade suction side-casing comer. Strong vorticity (negative from 
hub to mid span, and positive in the outer span) exists in the wake region. 

Further downstream at x/c, = 1.10 (see Fig. 37), the vortk ity structure remains nearly the same 
but at a lower level except in the wake close to the hub. The positive vorticity region associated with 
the leakage vortex at the blade suction side-casing comer becomes higher and is expanded at this axial 
location. The secondary vorticity seems to decay as flow progresses from x/q = 1 .01 to 1 . 10, but the 
leakage vorticity and the extent of this increases downstream. The vorticity in the wake-secondary 
flow interaction region decreases in strength and extent as the flow progresses downstream. 

4.5 Turbulence Intensity and Stresses 

Contours of total turbulence intensities, based on relative velocities, are shown in Fig. 38. 
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Turbulence intensity contours of components in the intrinsic co-ordinates (w s -, w„- , w r ) are shown in 
Fig. 19-21, and 39. These are normalized by the local total relative velocity. Very high values of 
turbulence intensities are measured in the tip vortex and secondary flow regions near the casing, and in 
the wakes. It is interesting to note that the radial components are much higher than the other two 
components, both in the wake and in the secondary and leakage flow regions. This is due to highly 
anisotropic nature of turbulence. The maximum turbulence levels in the casing secondary flow region 
is very high (25%), indicating highly turbulent nature of flow in the casing-suction surface corner. 
High turbulence intensities ( 15 - 18%) exist in nearly half of the outer passages. This is an indication 
that anisotropic turbulence models are essential to predict turbine endwall flows. 

In the casing half of the passage, there is a high turbulence level associated with the interaction 
between the tip leakage flow and the casing secondary flow. It influences most of the outer span. 
There is also a turbulence structure existing in the leakage-induced separated flow region. Along the 
top and bottom edges of the negative axial velocity region (Fig. 32, 37), the turbulence level is above 
30 - 35%. The remaining high turbulence distributions (around 20%) are in the two regions located 
about 15% pitch from the blade suction side and between H = 0.75 and 0.65 and between 0.40 and 
0.30 immersions. Through a re-examination of the secondary flows in Fig. 31, it is found that strong 
variations in the radial velocity occur exactly in these two regions. These two regions are also 
coincident with the severe wake distortion described in the previous section. Further downstream at 
x/c, = 1.10, the turbulence distribution remains the same, but in a lightly weaker form The structure 
associated with the separated flow is fully mixed out at this downstream location. 

The radial distribution of maximum, minimum and passage-averaged total turbulence intensities 
are plotted in Fig. 40 for x/c, =1.01 and 1.1. As indicated earlier (Figs. 31-35), no data was acquired in 
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the region H = 0.0 to 0.5 at x/Cr = 1.01 in the blade region due to inaccessibility by the laser. As 
mentioned earlier, the maximum intensities occur near the tip clearance vortex. This is caused by 
intense mixing between the secondary flow, the tip-leakage flow, and the blade and annulus wall 
boundary layers. In the casing secondary flow region, the intensities are high in most of the outer span. 
The minimum values occur in the inviscid core region, and its value ranges from 3 - 5%, based on the 
local velocity. 

The contour maps of the shear stresses are shown in Figures 41 through 44. Associated with 
the leakage vortex, there are a pair of high radial shear stress regions (w r - w n -) with opposite si gns, 
which increase in magnitude at x/c, = 1.10. It is interesting to note that the magnitude of r adial and 
streamwise stresses are of the same order of magnitude. This evidently indicates a non-isotropic nature 
of turbulence and invalidates the concepts based on isotropic eddy viscosity (constant value in all three 
directions). As indicated earlier, the curvature and rotation effects alter the structure of turbulence, 
with major influence in the radial direction. The r adial distribution of maximum , minimum, and average 
shear stresses (absolute levels) at x/Cr = 1.01, plotted in Fig. 45, show that the maximum value of r adial 
stresses are higher than the maximum value of streamwise stresses. This confirms the earlier 
conclusion that isotropic eddy viscosity models are not valid for computation of leakage and secondary 
flows. 

4.6 Rotor Wake Characteristics 

Due to the importance in understanding of rotor/stator interaction and the resulting vibration 
and noise, the features of the rotor blade wake are examined in this section. The total relative velocity 
distributions are used for characterization of the rotor wake features and for comparison with the 
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nozzle vane wake analysis performed in the same facility. As demonstrated in the previous section, 
there are very strong secondary flows and tip leakage flow downstream of the rotor. Their influence 
on the rotor wake is so powerful that, at some radial locations, the wake pattern is completely 
distorted. For these distorted wake patterns, it is found that the velocity defect and semi-width can 
not be determined unambiguously and these locations have been omitted from the plots. Near the hub 
(H=0.3 to 0.02), there is an extensive secondary flow region occupying almost the hub-half of the blade 
passage (see Fig. 31 and 37). This secondary flow increases the total relative velocity on the pressure 
side and decreases it on the suction side (e.g. Fig. 12 at H = 0. 10). The secondary flow sweeps the hub 
endwall flow from the pressure to the suction side, moving the low momentum fluid to the suction side 
and high momentum fluid to the pressure side. The degree of velocity reduction is so severe at H = 0.3 
and H = 0.2 immersions that the wake features (such as velocity defect and semi-width) are no longer 
identifiable. 

Similar phenomena are observed for the secondary flow region in the casing-half of the blade 
passage. For example, the wake pattern at H = 0.70 is so distorted (as distorted as at H = 0.30 radial 
location shown in Figs. 32 and 33) that no wake features can be identified. Very close to the casing, 
however, the influence of the tip clearance flow has a dominant effect (H = 0.99 through 0.93, Figs. 9- 
12, 32, 33). But, unlike the secondary flows, this influence is confined to the center of the wake; a 
much smaller influence region than that of the secondary flow. 

The velocity defect Wc/W^* and the semi-width (L/S) of the wake are plotted in Fig. 46. The 
main feature here is the relatively large increase in the velocity defect in the casing region due to the 
presence of the strong tip clearance vortex. The wake decay is much more rapid in the midspan region 
due to the thin, well defined wake present right downstream of the trailing edge of the blade. The 
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wake defect decreases dramatically (by as much as 50%) at the 10% downstream location. In the 
secondary flow region (H — 0.8) the wakes are not well defined due to the presence of strong passage 
vortex. The wake defect remains relatively constant at the 10% downstream location over most of the 
blade passage, due to the fact that most of the radial mixing in the wake has already occurred. As for 
the semi-width of the wake (Fig. 46b), there is an apparent increase in the width r adially outward. 
There are two local peaks at H = 0.75 and H = 0.4 as a result of the secondary flows, which tend to 
reduce the total relative velocity on the blade suction side and, thus, to widen the wake. There is an 
additional peak very close to the casing, which is caused by tie tip leakage flow. The secondary 
flow/wake interaction results in fester decay of the wakes in these regions. The semi-wake width at 
both locations is very s imilar 

The decay characteristics of the rotor wake velocity defect is plotted in Fig. 47 along with the 
decay features of the nozzle wake of the same stage (Zaccana and Lakshminarayana, 1995). In 
addition, some computational results and the experimental data for similar turning (100-110 degrees) 
cascade are also shown. The nozzle wake decays fester than the rotor wake. This is probably caused 
by the potential interaction (unsteady) between the nozzle and the downstream rotor. The rotor wake 
decay near the trailing edge (up to 10% of chord in the streamwise direction) is very rapid. Large 
radial velocities associated with the blade rotation tend to decay the wakes fester than the stationary 
blades as observed by Ravindranath and Lakshminarayana (1982). But the decay rate slows down 
downstream of the trailing edge. The decay of the rotor wake is in good agreement with the 
correlations of Raj and Lakshminarayana (1973). 

5. CONCLUDING REMARKS 

A 3-D LDV system is used to measure the unsteady three- dimensional flow field downstream 


22 



of a turbine rotor stage. Extreme caution is exercised to ensure reliable measurement of the radial 
component of velocity. Through extensive testing, the optimum angles for blue, violet , and green 
beams are established. Also, the window size and the minimum data number per window for the LDV 
data reduction process are tested, and suitable parameters are selected to avoid spurious fluctuations. 

The rotor wake pattern is much more disturbed by the three-dimensional flow field (secondary 
flows and tip leakage flow) than that observed in the nozzle wake. At some r adial locations, the 
distortion is so severe that the wake features (such as the velocity defect and the semi-width of the 
wake) can not be identified unambiguously. The transport of the nozzle wake and secondary flow 
through the rotor influences its flow field by increasing the suction surface boundary layer thickness 
and turbulence intensity, especially in the tip region. Unlike the nozzle wake, there is an increase in the 
rotor wake width in the radial direction due to rotation effects. From hub to mid-span, the secondary 
flow and rotation effects are cumulative in inducing radial outward flow in the wake. From mid-span 
to casing, the secondary flow effects dominate in inducing r adial inward velocity in the wake. 

The global features of the three-dimensional flow field downstream of the rotor stage are 
examined in terms of velocity, vorticity, turbulence, and shear stress distributions relative to the rotor 
blade. The hub secondary flow is found to be very strong, well organized, and occupying the entire 
blade passage. The casing-half secondary flow presents a much more fragmented pattern due to the 
interaction with the leakage flow. The secondary flow - rotor wake interaction region has low radial 
velocity, high axial vorticity, and high turbulence intensities (25 - 30%). 

The secondary flow and relative motion between the blade and the casing wall seem to reduce 
leakage flow, and the leakage vortex is confined to the comer of the suction surface - casing regions. 
Strong secondary flow, which tends to confine the tip vortex, reduces the diffusion and spreading of 
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the tip vortex downstream. The tip clearance vortex is confined to about 5% of span near the tip. 
There is a major difference in the behavior of the leakage vortex between those observed in cascades 
and in the present rotor. 

The leakage flow/vortex region is characterized by separated flow, high turbulence intensity 
and flow underturning. The secondary flow seem to overshadow all other effects from midspan to 
H = 0.95. The overall vorticity field due to secondary flow (both hub and tip) are much more 
dominant compared to the leakage vortex, which is confined to a very small region. 

The radial turbulence intensity was found to be substantially higher than other components in 
the secondary and leakage flows and wake regions due to the anisotropy of the turbulent flow field. 
Over most of the rotor blade passages the shear stress correlation terms are found to be very low 
except in the tip vortex and casing passage vortex regions. The maximum absolute values of wvvvv 
are higher than the maximum absolute values of w,-w r ' and w„- w s - indicating the need to use an 
isotropic eddy viscosity model in the computation of the rotor flow field. In the free stream of the 
rotor, a substantial increase in the turbulence level (3~8%) was observed compared to the nozzle free 
stream levels (1.5 ~ 2%), and this is caused by the interaction of the nozzle wake and the rotor flow 
field. 
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Figure 1. Schematic of the test section with locations of previous and current measurement 
locations. 


27 





(In meter) 


Figure 2. The grid of data points of the LDV measurement. The spatial resolution is increased in 
the tip region to study the leakage flow. 
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Figure 3. The blade periodic, blade aperiodic, revolution periodic, and revolution aperiodic 
components at the x/c,=1.01 (1% chord from trailing edge) and H=0.93 position. 
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Figure 5. The axial (a) velocity and (b) turbulence at 1% of chord downstream of the rotor and 
H=0.60 for the 50 and 100 data windows per blade passage. The axial velocity is chosen for the 
maximum variation, and the turbulence is almost identical for the two window sizes. 
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Figure 6. The variation of (a) data number per window and (b) radial velocity distribution with 
the radial location at 10% of chord downstream of the rotor. 
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Figure 7. The axial turbulence intensity distributions for 5 individual blade passages (based on 
ensamble averaged velocity) and for the averaged blade passaged. 
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Figure 8. Total relative velocity 1% downstream of the rotor at H=0.6. The vertical bars indicate 
the total unresolved unsteadiness at the measurement location. 
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Figure 9. Axial velocity distributions at 18 radial stations extending from H=0.02 through 0.99 at 
x/Cr=1.01. (Note that the blank regions are due to blade blockage of the laser beams at the 
measurement locations.) 


35 















Figure 10. Radial velocity distributions downstream of rotor (x/c,=1.01). Radial outward direction 
is positive: 
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Figure 11. (Continued) 












Figure 12. Total relative velocity distributions at x/c,=1.01. (Note that the blank regions are due 
to blade blockage to the LDV system in the measurement locations.) 













Figure 12. (Continued.) 



































Figure 14. Axial velocity distributions at 18 stations from H=0.02 through 0.99 at x/c,=1.10. 
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Figure 15. Radial velocity distributions downstream of rotor passage, H=0.02 through 0.99, at 
x/c,=1.10. 
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Figure 15. (Continued) 







































Figure 17. (Continued.) 



































































































































































Figure 28. (Cont.) 
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Figure 31. The secondary velocity vector (V^) maps at (a) 1% (x/c^l.01) and (b) 10% 
(x/Cr=1.10) of chord downstream of the rotor. 
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Figure 36. The comparison of the measured passage-averaged velocity at x/Cr=l . 10 with the 
design value. 
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(b) 

Figure 37. Measured Axial Vorticity (® X /Q) at (a) x/c^l.Ol and (b) x/c^l.10 downstream of the 
rotor. 
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Figure 39. Contours of measured turbulence intensities in the s,n,r coordinate system at x/Cr-1.01 
and 1.10. 
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Figure 45. Maximum absolute values of the shear stress distribution downstream of the rotor 
blade at (a) x/c,=1.01 and (b) x/Cr=1.10. 
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